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Abstract 

Multivariate versions of the Kronecker theorem in the continuous multivariate setting has recently 
been published. These theorems characterize the symbols that give rise to finite rank multidimen¬ 
sional Hankel and Toeplitz type operators defined on general domains. In this paper we study how 
the additional assumption of positive semidefinite affects the characterization of the corresponding 
symbols, which we refer to as Caratheodory-Fejer type theorems. We show that these theorems 
become particularly transparent in the continuous setting, by providing elegant if and only if state¬ 
ments connecting the rank with sums of exponential functions. We also discuss how these objects 
can be discretized, giving rise to an interesting class of structured matrices that inherit these desir¬ 
able properties from their continuous analogs. We describe how the continuous Kronecker theorem 
also applies to these structured matrices, given sufficient sampling. We also provide a new proof for 
the Caratheodory-Fejer theorem for block Toeplitz matrices, based on tools from tensor algebra. 
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1. Introduction 

The connection between low-rank Hankel and Toeplitz operators and matrices, and properties 
of the functions that generate them play a crucial role for instance in frequency estimation disni 
system identification dl [m US ED and approximation theory H El 1 [S IS dS EH- 
The reason for this is that there is a connection between the rank of the operators, and the fact 
that the functions that generate these operators and matrices, respectively, “generically” are sums 
of exponentials, where the number of terms is connected to the rank of the operators and matrices 
(Kronecker’s theorem). Moreover, adding the condition of positive semidefinite imposes further 
restrictions on the sums of exponentials (Caratheodory-Fejer’s theorem), a result which underlies 
e.g. Pisarenko’s famous method for frequency estimation [39]. 

We provide corresponding theorems in the multidimensional setting. In contrast to the one 
dimensional situation, the multidimensional framework provides substantial flexibility in how to 
define these operators. Whereas most previous research on multidimensional Hankel and Toeplitz 
type operators considers “symbols” or “generating sequences” / that are defined on product domains, 
we here consider a framework where / is defined on an open connected and bounded domain Q 
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in (satisfying mild assumptions). Besides providing beautiful new theorems, it is our hope 
that the new results in this paper will pave the way for applications in multidimensional frequency 
estimation/approximation/compression, in analogy with the use of Toeplitz and Hankel matrices 
in the one dimensional setting. For this reason, we present results both in the continuous and 
discretized setting, and discuss how they influence each other. 

To present the key ideas, we here focus mainly on the continuous theory since it is more trans¬ 
parent. Let us denote the multidimensional Hankel type operators by F^ and their Toeplitz coun¬ 
terparts by 0/ (see Section 2.2 for appropriate definitions). These operators were introduced in [3] 
where it is shown that if F/ or 0/ has rank K < oo, then / is necessarily an exponential polynomial; 


fix) = '^Pj(,x)e^x^ 

t=i 


( 1 . 1 ) 


where J < K, pj are polynomials in x = (xi,..., Xd), Cj ^ Q ■ x denotes the standard scalar 

product 

d 

Ck ‘ ^ ^ ^ 

m—1 


Conversely, any such exponential polynomial gives rise to finite rank Tj and 0 f respectively, and 
there is a method to determine the rank given the “symbol” (l.ll. Most notably, the rank equals 
K f is of the form 

K 


fix) = 




c^e 


Ck‘X 


( 1 . 2 ) 


where S C (assuming that there is no cancelation in (0)- 

In this paper, we study what happens if one adds the condition that Tf or Qf be positive 
semi-definite (PSD). We prove that 0y then has rank K if and only if / is of the form 


K 

/(x)=^Cfce*«-" (1.3) 

k^l 

where Cfc > 0 and G (Theorem |7.l| , which in a certain sense is an extension of Caratheodory- 
Fejer’s theorem on PSD Toeplitz matrices. Correspondingly, Tf is PSD and has rank K if and only 
if / is of the form 

K 

fix) =^Cke^^'^ (1-4) 

k^l 

where again > 0 and G (Theorem |8.l| . Similar results for Hankel matrices date back to 
work of Fischer m- 

The only of the above results that has a simple counterpart in the finite dimensional discretized 
multivariable setting is the Caratheodory-Fejer’s theorem, which has been observed previously in 
|49] (concerning block Toeplitz matrices). In this paper we provide a general result on tensor 
products, which can be used to “lift” structure results in one-dimension to the multi-dimensional 
setting. We use this to give an alternative proof of the discrete Caratheodory-Fejer theorem, which 
subsequently is used to prove the continuous counterpart. 
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Figure 1: a) “Generating sequence” defined on a disc b) The matrix realization of the corresponding “general 
domain Hankel operator” (see Section [4.1| for further details). 


Fischer’s theorem on the other hand has no neat version in the multivariable finite dimensional 
setting, but has been generalized to so called small Hankel operators on in [ID], a paper 

which also contains a result analog to (1.4|. 

However, the product domain setting is rather restrictive and not always a natural choice. 
Whereas one-dimensional “symbols” necessarily are defined on an interval, there is an abundance of 
possible regions to define their corresponding multidimensional cousins. Despite this, the majority 
of multivariate treatments of these issues are set either directly in a block-Toeplitz/Hankel setting, 
or rely on tensor products. In both cases the corresponding domain of definition D of the symbol is 
a square (or multi-cube), but for concrete applications to multidimensional frequency estimation, 
the available data need not be naturally defined on such a domain. In radially symmetric problems, 
a circle may be more suitable or, for boundary problems, a triangle might be more appropriate. 

Concerning analogs of the above results for the discretized counterparts of 0/ and F j, we show in 
this paper how to construct “structured matrices” that approximate their continuous counterparts, 
and hence can be expected to inherit these desirable properties, given sufficient sampling rate. We 
give simple conditions on the regularity of / and il needed for this to be successful. This gives 
rise to an interesting class of structured matrices, which we call “general domain Hankel/Toeplitz 
matrices”. As an example, in Figure [l] we have a “generating sequence” / on a discretized disc, 
together with a plot of its general domain Hankel matrix. 

The paper is organized as follows. In the next section we review the theory and at the same time 
introduce the operators we will be working with in the continuous setting (Section |2.2[ ) . The short 
Section provides a tool from tensor algebra, and also introduce useful notation for the discrete 
setting. Section discuss how to discretize the T/’s and 0/’s, and we discuss particular cases such 
as block Toeplitz and Hankel matrices. In Section]^ we prove the Caratheodory-Fejer theorem in 
the discrete (block) setting. Section]^ shows that the discrete operators approximate the continuous 
counterparts, given sufficient sampling rate, and we discuss Kronecker’s theorem. Sections and 
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considers structure results under the PSD condition, first for 0/’s and then for P/’s. In the last 
section, we extend the above results to the corresponding operators on unbounded domains. 


2. Review of the field 


Suppose that a Hankel matrix iJ or a Toeplitz matrix T of size N x N is taken from samples of 
a “discretized exponential polynomial” 

,7 

( 2 . 1 ) 

i=i 

(where Xj G C and n = 1,..., 2N — 1) of total degree 


K = ^(degpj + 1) 
i=i 


( 2 . 2 ) 


strictly less than N. Based on the fundamental theorem of algebra, one can show that the rank of 
either H or T equals K, and that the polynomials pj and the Xj’s are unique. This observation is 
sometimes used hand-wavingly in the converse direction, which is not true. However, in terms of 
applications this doesn’t matter because of the following stronger statement: If T or iP has rank 
K < N then its generating sequence is “generically” of the form 

K 

Y^cuXl- (2.3) 

/c=l 

In the multidimensional setting, the corresponding statement is false. We refer to [3] for a longer 
discussion of these issues, especially Section 8. See also goilis]. As an example of an exceptional 
case in the one-variable situation, consider the Hankel matrix 


/ 1 0 0 0 0 \ 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

\ 0 0 0 0 1 / 


(2.4) 


Clearly, the rank is 2 but the generating sequence (1, 0, 0,0,0,0, 0, 0,1) is neither of the form (2.1 \ 
nor (2.31. The book [^, which has two sections devoted entirely to the topic of the rank of 
finite Toeplitz and Hankel matrices, give a number of exact theorems relating the rank with the 
“characteristic” of the corresponding matrix, which is a set of numbers related to when determinants 
of certain submatrices vanish. Another viewpoint has been taken by B. Mourrain et. al dniiii 
l32l [33] , in which, loosely speaking, these matrices are analyzed using projective algebraic geometry 
and the 1 to the bottom right corresponds to the point oo. 

In either case, the complexity of the theory does not reflect the relatively simple interaction 
between rank and exponential sums, as indicated in the introduction. There are however a few 
exceptions in the discrete setting. Kronecker’s theorem says that for a Hankel operator (i.e. an 
infinite Hankel matrix acting on ^^(N)), the rank is K if and only if the symbol is of the desired 
form (0° defined as 1), with the restriction that |Ajj < 1 if one is only interested in bounded 
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operators, see e.g. [151 [?7l 1^157], Also, it is finite rank and PSD if and only if the symbol is of the 
form ( |2.3| ) with Cfe > 0 and A/c S (—1,1), a result which also has been extended to the multivariable 
(tensor product) setting |4fl| . In contrast, there are no finite rank bounded Toeplitz operators (on 
t!^(N)). If boundedness is not an issue, then a version of Kroneckers theorem holds for Toeplitz 
operators as well |18| . 

Adding the PSD condition for a Toeplitz matrix yields a simple result which is valid (without 
exceptions) for finite matrices. This is the essence of what usually is called the Caratheodory-Fejer 
theorem. The result was used by Pisarenko |39| to construct an algorithm for frequency estimation. 
Since then, this approach has rendered a lot of related algorithms, for instance the MUSIC method 
|43| . We reproduce the statement here for the convenience of the reader. For a proof see e.g. 
Theorem 12 in |5] or Section 4 in |21|. Other relevant references include HJdS]. 

Theorem 2.1. Let T be a finite TV + 1 x + 1 Toeplitz matrix with generating sequence {fn)n--N- 
Then T is PSD and Rank T = K < N if and only if f is of the form 

K 

f{n)=Y,cuK (2.5) 


where > 0 and the Xk ’s are distinct and satisfy |Afc| = 1. 


The corresponding situation for Hankel matrices H is not as clean, since (2.4) is PSD and has 
rank 2, but do not fit with the model (2.5) for > 0 and real A^’s. Results of this type seems to go 


back to Fischer m, and we will henceforth refer to statements relating the rank of PSD Hankel- 
type operators to the structure of their generating sequence/symbol, as “Fischer-type theorems” 
(see e.g. Theorem 5 |1] or |5T]). Corresponding results in the full rank case is found e.g. in |4()| . 

We end this subsection with a few remarks on the practical use of Theorem |2.1| For a finitely 
sampled signal, the autocorrelation matrix can be estimated by H*H where H is & (not necessarily 
square) Hankel matrix generated by the signal. This matrix will obviously be PSD, but in general 


it will not be Toeplitz. However, under the assumption that the Afc’s in (2.5) are well separated. 


the contribution from the scalar products of the different terms will be small and might therefore 
be ignored. Under these assumptions on the data, the matrix H*H is PSD and approximately 
Toeplitz, which motivates the use of the Caratheodory-Fejer theorem as a means to retrieve the 
Afc’s. 


2.1. Finite interval convolution and correlation operators 

The theory in the continuous case is much “cleaner” than in the discrete case. In this section we 
introduce the integral operator counterpart of Toeplitz and Hankel matrices, and discuss Kronecker’s 
theorem in this setting. 

Given a function on the interval [—2,2], we define the finite interval convolution operator : 

'df{g){x) = J f{x-y)g{y) dy. (2.6) 

Replacing x — y hy x + y we obtain the finite interval correlation operators T f. These operators 
also go by the names Toeplitz and Hankel operators on the Paley-Wiener space. It is easy to see 
that if we discretize these operators, i.e. replace integrals by finite sums, then we get Toeplitz and 
Hankel matrices, respectively. More on this in Section |4T| 
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Kronecker’s theorem (as formulated by Rochberg lU) then states that Rank 0y = K (and 
Rank F/ = K) if and only if / is of the form 


f{x) = (2.7) 

i=i 

where pj are polynomials and Q € C. Moreover, the rank of 0/ (or F/) equals the total degree 

j 

K = ^(degpj + 1). (2.8) 

f=i 

However, functions of the form 


K 


: Ck € ^ 


(2.9) 


are known to be dense in the set of all symbols giving rise to rank K finite interval convolution 


operators. Hence, the general form (2.71 is hiding the following simpler statement, which often is 


of practical importance. 0/ generically has rank K if and only if f is a sum of K exponential 
funetions. As already noted, this is false in several variables, see [3]. The polynomial factors appear 


in the limit if two frequencies in (2.9) approach each other and interfere destructively, e.g. 


X = lim 
e—>0+ 


-1 


( 2 . 10 ) 


This can heuristically explain why these factors do not appear when adding the PSD condition. 


since the functions on the right of (2.10) give rise to one large positive and one large negative 
eigenvalue. 


2.2. The multidimensional eontinuous setting: TCO’s 

Given any (square integrable) function / on an open connected and bounded set 12 in . 


d>l, 

the natural counterpart to the operator (2.6) is the (General Domain) Truncated Gonvolution 
Operator (TGO) 0/ : L^(T) —)■ T^(S) defined by 


%(5)(®) = / f{x-y) 9 {y) dy, x£ 
Jt 


where ^ and T are connected open bounded sets such that 

= E —'T = {x — y: x G E, y G T}. 


( 2 . 11 ) 


( 2 . 12 ) 


In [3] such TGO operators are studied, and their finite rank structure is completely characterized. 
It is easy to see that 0y has rank K whenever / has the form 


K 


f(x) a; G 12, 


(2.13) 


fc = l 
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where the (^i,... ,(k S C"^. 

The reverse direction is however not as neat as in the one-dimensional case. It is true that 
the rank is finite only if / is an exponential polynomial (Theorem 4.4 in [5), but there is no 
counterpart to the simple formula (2.8l. However, Section 5 (in 0 ) gives a complete description 
of how to determine the rank given the symbol / explicitly, Section 7 gives results on the generic 
rank based on the degree of the polynomials that appear in /, and we also provide lower bounds, 
and Section 8 investigates the fact that polynomial coefficients seem to appear more frequently in 
the multidimensional setting. Section 9 contains an investigation related to boundedness of these 
operators for the case of unbounded domains, which we will treat briefly in Section]^ of the present 
paper. 

If we instead set H = S -f T then we may define the Truncated Correlation Operator 


r/(5)(a^) = / f{x + y)g{y) dy, xG 


(2.14) 


This is the continuous analogue of finite Hankel (block) matrices. As in the finite dimensional case, 
there is no real difference between Tf and Of regarding the finite rank structure. In fact, one 
turns into the other under composition with the “trivial” operator c{f){x) = f{—x), and thus all 
statements concerning the rank of one can easily be transferred to the other. 


2.3. Other multidimensional versions 

The usual multidimensional framework is that of block-Hankel and block-Toeplitz matrices, 
tensor products, or so called “small Hankel operators on £^(N‘^). In all cases, the generating sequence 
/ is forced to live on a product domain. For example, in |48| they consider the generating sequences 


of the form (1.21 (where cc is on a discrete grid) and give conditions on the size of the block Hankel 
matrices under which the rank is K, and in it is observed that the natural counterpart of the 
Caratheodory-Fejer theorem can be lifted by induction to the block Toeplitz setting. For the full 
rank case, factorizations of these kinds of operators have been investigated in |20L 145] . Extensions 
to multi-linear algebra are addressed for instance in [5411551155] . 

Concerning “small Hankel operators”, in addition to [40] we wish to mention [25] where a formula 
for actually determining the rank appears, although this is based on reduction over the dimension 
and hence not suitable for non-product domains. 

There is some heuristic overlap between |3| and |221 [55] . In they consider block Hankel 
matrices with polynomial symbols, and obtain results concerning their rank (Theorem 4.6) that 
overlap with Propositions 5.3, Theorem 7.4 and Proposition 7.7 of |5| for the 2d case. Proposition 
7 in |25| is an extension to 2d of Kronecker’s theorem for infinite block Hankel matrices (not 
truncated), which can be compared with Theorem 4.4 in |5|. 

In the discrete setting, the work of B. Mourrain et al. considers a general domain context, 
and what they call “quasi Toeplitz/Hankel matrices” correspond to what here is called ’’general 
domain Toeplitz/Hankel matrices” (we stick to this term since we feel it is more informative for 
the purposes considered here). See e.g. Section 3.5 in [55], where such matrices are considered for 
solving systems of polynomial equations. In m, discrete multidimensional Hankel operators (not 
truncated) are studied, and Theorem 5.7 is a description of the rank of such an operator in terms 
of decompositions of related ideals. Combined with Theorem 7.34 of HZ], this result also implies 
that the symbol must be of the form (See also Section 3.2 of [3^ , where similar results are 

presented.) These results can be thought of as a finite dimensional analogue (for product domains) 
of Theorem 1.2 and Proposition 1.4 in |5]. Theorem 5.9 gives another condition on certain ideals in 
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order for the generating sequence to be of the simpler type, i.e. the counterpart of (1.21 instead of 
In Section 6 of the same paper they give conditions for when these results apply also to the 
truncated setting, based on rank preserving extension theorems. Similar results in the one-variable 
setting is found in Section 3 of |18| . 

Finally, we remark that the results in this paper concerning finite rank PSD Hankel operators 
partially overlap heuristically with results of m and those found in Section 4 in [3^, where the 
formula (2.31 is found in the (non-truncated) discrete environment. In the latter reference they 
subsequently provide conditions under which this applies to the truncated setting. 

With these remarks we end the review and begin to present the new results of this paper. For 
the sake of introducing useful notation, it is convenient to start with the result on tensor products, 
which will be used to “lift” the one-dimensional Caratheodory-Fejer theorem to the multidimensional 
discrete setting. 


3. A property of tensor products 

Let Ui,... ,Ud be finite dimensional linear subspaces of C”. Then 0j^iUj is a linear subspace 
of and the latter can be identified with the set of C-valued functions on {1,..., 

Given / G and x G {!,... we will write f(x) for the corresponding value. For fixed 

X = (xi, ..., Xd-i) S {1, ..., we define 

flM = . . . ,Xd),f 2 (x) = f(xi,-,X2,-- ■ ,Xd), h{x) =/(xi,X2,-,X3,-- - ,Xd), 

i.e. the vectors obtained from / by freezing all but one variable. We refer to these vectors as 
“probes” of /. If f G then it is easy to see that all probes fj of / will be elements of Uj, 

j = 1,... ,d. The following theorem states that the converse is also true. 

Theorem 3.1. If all possible probes fj{x) of a given f G lie in Uj, then f G 0j^iUj. 

Proof. First consider the case d = 2. Let V C Gj^iC" consist of all / with the property stated in 
the theorem. This is obviously linear and ?7i G C V. If we do not have equality, we can pick an 
f in V which is orthogonal to [/i G J72- At least one probe /i(-) must be a non-zero element ui of 
Ui. Given any U 2 G U 2 consider 

n n 

{ui G U 2 , f) = (u2,^Ml(i)/2(f)). (3.1) 

i=i 

From the middle representation and the choice of iti, we see that at least one value of the vector 
^i(*)/ 2 (*) is non-zero. Moreover this is a linear combination of probes / 2 (f)j and hence an 
element of U 2 . But then we can pick U 2 £ U 2 such that the scalar product (13 is non-zero, which 
is a contradiction to the choice of /. The theorem is thus proved in the case d = 2. 

The general case now easily follows by induction on the dimension, noting that can be 

identified with G {®'j= 2 ^'^) and that Gj^it/j under this identification turns into Ui G {®‘j^ 2 ^ 3 )- 

□ 


4. General domain Toeplitz and Hankel operators and matrices 

The operators in the title arise as discretizations of general domain truncated convolution/correlation 
operators. These become “summing operators”, which can be represented as matrices in many ways. 





Matrix realization of 0/ 
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Figure 2: Left: Domains S, T, and f2 = H — T with lexicographical order. Right: Illustration of where the 
numbered points in fl show up in the corresponding matrix realization of © j. 


4 ..I. Discretization 

For simplicity of notation, we here discretize using an integer grid, since grids with other sam¬ 
pling lengths (these are considered in Section 6.1) can be obtained by first dilating the respective 
domains. We will throughout the paper use bold symbols for discrete objects, and normal font for 
their continuous analogues. Set 

T = {x e : X e T}, 


make analogous definition for S/H and define SI = Y — H. We let ©j denote what we call a “general 
domain Toeplitz operator” 


®f{ 9 ){x)=^ f{x-y)g{y), x e H, (4.1) 

where g is an arbitrary function on T. We may of course represent 5 as a vector, by ordering the 
entries in some (non-unique) way. More precisely, by picking any bijection 


o,:{l,...,|Y|}^ Y 


(4.2) 


we can identify g with the vector g given by 

(5i)i=i = gioyU))- 

Letting Ox be an analogous bijection for H, it is clear that 0/ can be represented as a matrix, 
where the (i, j)’th element is f{ox{i) — Oy{j)). Such matrices will be called “general domain Toeplitz 
matrices”, see Figurej^for a small scale example. We make analogous definitions for F f and denote 
the corresponding discrete operator by Fy. We refer to this as a “general domain Hankel operator” 
and its matrix realization as “general domain Hankel matrix”. An example of such is shown in 
Figure [2 
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^.2. Block Toeplitz and Hankel matrices 


If we let H and T be multi-cubes and the ordering bijections be the lexicographical order, 
then the matrix realization ©/ of (4.1) becomes a block Toeplitz matrix. These are thus special 
cases of the more general operators considered above. Similarly, block Hankel matrices arise when 
representing Tin the same way. 

For demonstration we consider H = Y = { — 1,0,1}^ so SI = {—2,..., 2}^. The lexicographical 
order then orders { — 1 , 0 , 1 }^ as 


( 1 , 1 , 1 ), ( 1 , 1 , 0 ), ( 1 , 1 ,- 1 ), ( 1 , 0 , 1 ), ( 1 , 0 , 0 ),...,(- 1 ,- 1 ,-!). 


The matrix-realization T of a multidimensional Toeplitz operator ©f then becomes 


/^/3(0 

0) 

^/3(0 -1) 

'^/3(0.- 

-2) 

^/3(-1.0) 


^/3(- 

1-2) 

^/3(- 

2,0) 

'^/3(- 

2,-1) 

^/3(- 

2,-2) 

\ 

^/3(0 

1) 

^/3(0.0) 

'^/3(0.- 

-1) 

^/3(-l,l) 

^/3(-1.0) 

^/3(- 

1-1) 

^/3(- 

2.1) 

'^/3(- 

2,0) 


2.-1) 


^/3(0 

2) 

^/3(0.1) 

^/3(0.0) 

^/3(-1.2) 

^/3(-l.l) 

^/3(- 

1.0) 

Tf3(- 

2,2) 

'^/3(- 

2.1) 

Tf3(- 

2,0) 


^/3(1 

0) 

^/3(1 -1) 

^/3(1.- 

-2) 

^/3(0.0) 

3(0-1) 

^/3(0 

-2) 

Tf3(- 

1.0) 

'^/3(- 

l.-l) 

Tf3(- 

1.-2) 


^/3(1 

1) 

^/3(1,0) 

'^/3(1.- 

-1) 

^/3(0,1) 

3(0,0) 

^/3(0 

-1) 

Tf3(- 

1.1) 

'^/3(- 

1.0) 

Tf3(- 

l.-l) 


^/3(1 

2) 

^/3(1,1) 

^/3(1.0) 

^/3(0.2) 

^/3(0.1) 

^/3(0 

0) 

Tf3(- 

1.2) 

'^/3(- 

1.1) 

Tf3(- 

1.0) 


^/3(2 

0) 

^/3(2 -1) 

^/3(2.- 

-2) 

^/3(1,0) 

^/3(1.-1) 

^/3(1 

-2) 

^/3(0 

.0) 

'^/3(0 

-1) 

^/3(0 

-2) 


^/3(2 

1) 

^/3(2.0) 

^/3(2,- 

-1) 

^/3(1,1) 

^/3(1.0) 

^/3(1 

-1) 

^/3(0 

.1) 

'^/3(0 

0) 

^/3(0 

-1) 


V/3(2 

2) 

^/3(2.1) 

^/3(2.0) 

^/3(1.2) 

^/3(1.1) 

^/3(1 

0) 

^/3(0 

.2) 

'^/3(0 

1) 

^/3(0 

0) 

/ 


where e.g. 


/ /( 0 , 0 , 0 ) 

T’/3(o,o)= /(0,0,1) 

V /(0,0,2) 


/(0,0,-l) /(0,0,-2) \ 

/(0,0,0) /(0,0,-l) 

/( 0 , 0 , 1 ) /( 0 , 0 , 0 ) ) 


Note that this matrix has a Toeplitz structure on 3 levels, since each 3 x 3-block of the large matrix 
above is Toeplitz, and these blocks themselves form a 3 x 3 Toeplitz structure. 


4 ..3. Exponential sums 

We pause the general development to note some standard facts that will be needed in what 
follows. Fix iV € N, and for j = 1,..., d let be a set of at most 2N numbers in C. Set 
$ = $ix,... X (hd. 

Proposition 4.1. The set {e^'* : C G $} zs linearly independent as functions on {—N,... ,N}‘^. 

Proof. If d = 1 the result is standard, see e.g. Proposition 1.1 in |TH] or [12 Sec. 3.3]. For d > 1, 
the function e‘''^ = is a tensor. The desired conclusion now follows from the d = 1 

case and standard tensor product theory. □ 

We now set Y = H = {—N ,..., N}‘^, and let = {—2N ,..., 2^^}^^ in accordance with subsec¬ 
tion |4.1| Consider functions / on ft given by 


K 

/(®) = 

k^l 


(4.3) 


We say that the representation (4.3) is reduced if all Cfc’s are distinct and the corresponding coeffi¬ 


cients Ck are non-zero. 
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Proposition 4.2. Let $ be as before. Let the function f on {—2N,... ,2N}'^ be of the reduced 
form (4.3) where each fk is an element o/$. Then 


Rank 0/ = Rank P/ = K. 

Proof. Pick a fixed ( and consider f{x) = then 




®/(5)(a^) = 

yer 

which has rank 1. F or a general / of the form ( |4.3[ ) the rank will thus be less than or equal to K. 
But Proposition 4.1 implies that the set {e‘’'' ®}^7^is linearly independent as functions on H. Thus 
the rank will be precisely K, as desired. □ 


We end this section with a technical observation concerning 1 variable. 

Proposition 4.3. Let f be a vector of length m > n + 1 and K < n. Let be fixed and 

suppose that each sub-vector of f with length n + 1 can be written of the form (4.3), then f can be 
written in this form as well. 

Proof. Consider two adjacent sub-vectors with overlap of length n. On this overlap the representa¬ 
tion (4.3) is unique, due to Proposition 4.1 The result now easily follows. □ 


5. The multidimensional discrete Caratheodory-Fejer theorem 

Throughout this section, let T, H and be as in Sections |4.2| and [4.3[ i.e. multi-cubes centered 
at 0. The following theorem was first observed in [JS], but using a completely different proof. 

Theorem 5.1. Given f on {—2N,... ,2N}‘^, suppose that &f is PSD and has rank K where 
K < 27V. Then f can be written as 

K 

f{x) = (5.1) 

k=l 

where Ck > 0 and are distinct and unique. Conversely, if f has this form then &f is PSD 

with rank K. 


Proof. First assume that ©y is PSD and has rank K. Let T be a block Toeplitz representation 
of ©/, as described in Section 4.2 
decomposition as 


Since the rank of T is K, we can write its singular value 


K 


T = ^ SkUkvl. 


(5.2) 


k=l 


where Sk > 0, Uk, Vk are column matrices and T denotes transpose. Recall from Sectionthat the 


Toeplitz matrix 7/^(0) is the 2N -fix 2N 1 sub-matrix on the diagonal of T, (and 0 S ^). 


This is clearly PSD so by the classical Caratheodory-Fejer theorem (Theorem 2.1), fd{0) can be 
represented by 


K 


a; e {-27V... 27V}. 


(5.3) 


k^l 
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We identify functions on {—N... N} with in the obvious way, and define Ui C by 

Ui = Span 


The analogous subspace of will be called 


Note that /d(0) e 


Uext 


by (5.3). Set 


^■1 = {Cl) ■ • ■ )Ck}- 

For 1 < m < (27V + 1)‘^“^, let Uk^m denote the sub-vector of Uk related to the mth column-block, 
(i.e. with subindex ranging between (m — l)(2fV-|-1) -I-1 and m(2N -1-1)), and define Vk,m similarly. 


Then (5.2) implies that 


K 


T, 


/d(0) 




■'k^m' 




But this means that each Uk^m is in Ui, since /d(0) G Uf^*. Fix y G {—2N,..., 2iV}‘^ Restricting 


(5.2) to a corresponding off-diagonal Toeplitz matrices in T gives, with appropriate choice of m and 


j, the representation 


K 




fd{y) 


= E 


SkUk,m'Vk^rn+j- 


k=l 


But this means that the columns of Tf^ (y) l ie in C/i. We conclude that each sub-vector of fd{y) of 
length 2N -|- 1 is in Ui. By Proposition 4.3 we conclude that each probe fd{y) is in [/“*. 


By choosing a different ordering and repeating the above argument, we conclude that for each I 
{1 < I < d), there is a corresponding subspace such that all possible probes /;(•) are in 
Let ^k G be an enumeration of all IC‘^ multi-frequencies $i x ... x The corresponding K‘^ 


exponential functions span By Theorem 3.1 we can thus write 


fix) = 


E 


Ck^ 




(5.4) 


However, by Proposition 4.2 precisely K of the coefficients Ck are non-zero. This is (|5.1[). The 


uniqueness of the multi-frequencies is immediate by Proposition 4.1 (applied with N := 2N). The 
linear independence of these functions also give that the coefficients are unique. To see that Ck is 
positive, (1 < fc < K), just pick a multi-sequence on H which is orthogonal to all other j yf k. 


Using the representation (4.1) it is easy to see that 


O<(0/(5),ff)=Cfe|(g,e*«'“-)|2 


(5.5) 


and the first statement is proved. 

For the converse, let / be of the form (5.1). Then ©y has rank K by Proposition 4.2 and the 
PSD property follows by the fact that 


K 


O<(0/(5),-7) = E^^K5,e*«'='")| = 


(5.6) 


fc=i 


in analogy with (5.5). 


□ 


It is possible to extend this result to more general domains as considered in Section |4.1[ How¬ 
ever, such extensions are connected with some technical conditions, which are not needed in the 
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continuous case. Moreover, in the next section we will show that the discretizations of Section [4T] 
capture the essence of their continuous counterparts, given sufficient sampling. For these reasons 
we satisfy with stating such extensions for the continuous case, see Section 

The above proof could also be modified to apply to block Hankel matrices, but since Fischer’s 
theorem is connected with preconditions to rule out exceptional cases, the result is not so neat. (It 
does however provide alternative proofs the results in |in],i.e. concerning small Hankel operators). 
We here present only the cleaner continuous version, see Section 


6. The multidimensional discrete Kronecker theorem 


If we want to imitate the proof of Theorem Eii in Kronecker’s setting, i.e. without the PSD 
assumption, then we have to replace (5.3) with (2.7). With suitable modifications, the whole 
argument goes through up until (5.4), where now the ^^’s can lie in C'^ and Cfc also can be monomials. 


However, the key step of reducing the ^-term) representation (5.4) to the (7f-term) representation 
(5.11, via Proposition 4.2 fails. Thus, the only conclusion we can draw is that / has a representation 


fix) = ^pj{x)e^ 
i=i 


a: G D, 


( 6 . 1 ) 


where J < K, but we have very little information on the amount of terms in each pj. This is a 
fundamental difference compared to before. In examples are presented of TCO’s generated by 
a single polynomial p, where Fp has rank K much lower than the amount of monomials needed to 
represent p. It is also not the case that these polynomials necessarily are the limit of functions of 


the form (2.13) (in a similar way as (2.10)), and hence we can not dismiss these polynomials as 


“exceptional”. To obtain similar examples in the finite dimensional setting considered here, one can 


just discretize the corresponding TCO’s found in |3] (as described in Section 4.1). 


Nevertheless, in the continuous setting (i.e. for operators of the form 0y and Fj, c.f. (2.11) and 


(2.14)) the correspondence between rank and the structure of / is resolved in |3]. In particular it is 
shown that (either of) these operators have finite rank if and only if / is an exponential polynomial, 
and that the rank equals iF if / is of the (reduced) form 


K 

f = 


( 6 . 2 ) 


We now show that these results apply also in the discrete setting, given that the sampling is 
sufficiently dense. For simplicity of notation, we only consider the case F f from now on, but include 
the corresponding results for 0/ in the main theorems. 

6.1. Discretization 

Let bounded open domains T, S be given, and let Z > 0 be a sampling length parameter. Set 

T; = G T : n G Z'^}, 


(c.f. (4.1)), make analogous definition for H/ and define fi; = T/ + H;. We denote the cardinality 


of Y; by |T;|, and we define ^^(T;) as the Hilbert space of all functions g on T; and norm 

Mp = E \9iy)\^- 

ye-Tt 
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We let Tfj : —>■ denote the summing operator 

^fA9){x)= fix + y)g{y), xgSi. 

yeTi 

When I is understood from the context, we will usually omit it from the notation to simplify the 
presentation. It clearly does not matter if / is defined on S + T or H/ + Tp and we use the same 
notation in both cases. We define 0 fj in the obvious analogous manner. Note that in Section 
and we worked with 0y, which with the new notation becomes the same as 0/.i. 

Proposition 6.1. There exists a constant C > 0 such that 

Proof. By the Cauchy-Schwartz inequality we clearly have 

|r'/,;(5')(a;)| < \\fsni\\p{nt)\\9\\e^{-rt) 

for each x G Si. If we let |Hi| denote the amount of elements in this set, it follows that 

l!r/,;|| < 

Since S is a bounded set, it is clear that is bounded by some constant, and hence the result 

follows. □ 

Theorem 6.2. Let f G L^(fl) be continuous. Then 

Rank T/ ; < Rank F/ 


(and Rank 0y ; < Rank 0/J. 

Proof. Given y G Ti and f < Z let denote the multi-cube with center y and sidelength t, i.e. 
C'y* = {y GM.'^ : \y — y\ao < Z/2}) where | • |oo denotes the supremum norm in Chose to such that 
v^to/2 < dist(T;,9T). For t <to we then have that the set {e^‘}ygXi defined by 
is orthonormal in L^(T). We make analogous definitions for H;. Clearly A{Ti) is in bijective 
correspondence with Span {e^‘}ygT, via the canonical map P*’*, i.e. P^'*{5y) = where 5y is the 
“Kronecker 5—function”. Let denote the corresponding map : A(Si) -G Lf[S). 

Now, clearly Rank < Rank Tf and 


1 

¥ 


(Q'’* TfP^ 


Lt 




1 

f2d 


\x-x\^<t/2 J\v-v\aa<t/2 


f{x + y) dy dx. 


If we denote this number by /*(a; + y), we see that ^Q*’**r fP^'^ = Fjt It follows that Rank F^t ^ < 
Rank F f. Since / is continuous, it is easy to see that lim(_,,o+ f*{x + y) = f{x + y), which implies 
that limj_,,o+ ft i = Tf,i, and the proof is complete. □ 
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6.2. From discrete to continuous 

Our next result says that for sufficiently small I, the inequality in Theorem |6.2| is actually an 
equality. This needs some preparation. Given y G Ti we now define Cy to be the multi-cube 
with center y and sidelength 1. Set T™* = {y G T/ : Cy C T}, i.e. the set of those y’s whose 
corresponding multicubes are not intersecting the boundary. Moreover, for each y G T;, set 

, ^ r , */yGTr‘ 

^ \ 0, else 

We now define P* : £^(T;) —>■ L^(T) via P^{5y) = e^. Note that this map is only a partial isometry, 
in fact, pP is the projection onto Span : y G T^*}, and p^pP is the projection in L^(T) 
onto the corresponding subspace. We make analogous definitions for H/, denoting the corresponding 
partial isometry by QK Set 

Ni=Ni{T) = |Tz\Tr‘|, 

i.e. Ni is the amount of multi-cubes Cy intersecting the boundary of T, and note that Ni = 
diniKer PK Since T is bounded and open, it is easy to see that |T;"*| is proportional to We 
will say that the boundary of a bounded domain T is well-behaved if 

lim 1‘^Ni = 0. (6.3) 

1 -^ 0 + 

In other words, dT is well behaved if the amount of multi-cubes Cy properly contained in T 
asymptotically outnumbers the amount that are not. 

Proposition 6.3. Let T be a bounded domain with Lipschitz boundary. Then dT is well behaved. 

Proof. By definition, for each point x G i9T one can find a local coordinate system such that dT 
locally is the graph of a Lipschitz function from some bounded domain in to K, see e.g. m 

or m, Sec. 4.2. It is not hard to see that each such patch of the boundary can be covered by a 
collection of balls of radius I, where the amount of such balls is bounded by some constant times 
Since i9T is compact, the same statement applies to the entire boundary. However, it is 
also easy to see that one ball of radius I can not intersect more than 3^^ multi-cubes of the type 
Cy, and henceforth TV; is bounded by some constant times l/l'^~^ as well. The desired statement 
follows immediately. □ 

We remark that all bounded convex domains have well behaved boundaries, since such domains 
have Lipschitz boundaries, (see e.g. m Sec. 6.3]). Also, note that the above proof yielded a faster 
decay of Nil‘^ than necessary, so most “natural” domains will have well-behaved boundaries. We 
are now ready for the main theorem of this section: 

Theorem 6.4. Let the boundaries of T and S be well behaved, and let f be a continuous function 
on cl{Ll). Then 

Tf = lim l^^Q^FfjP^*. (6.4) 

1 ^ 0 + 

(Also 0/ = lim/_,,o+ l^Q''®f^iP^*). 
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Proof. We first establish that p^p^* converge to the identity operator I in the S'OT-topology. Let 
g G L^(T) be arbitrary, pick any e > 0 and let 5 be a continuous function on cZ(T) with ||g —g|| < e. 
Then 

||5 - P^P^*9\\ < \\9 - ~9\\ + \\9 - P^P^*~9\\ + \\P^P^\~9 - 9)\\- 

Both the first and the last term are clearly < e, whereas it is easy to see that the limit of the middle 
term as I ^ 0^ equals 0, since g is continuous on cl{T) and the boundary is well-behaved. Since e 
was arbitrary we conclude that lim/_>g+ p’-p’- g = g, as desired. The corresponding fact for is 
of course then also true. 

Now, since Fj is compact by Corollary 2.4 in [5], it follows by the above result and standard 
operator theory that 


•I vl* 


Tf = lim TfP'P 

1 ^ 0 + 


and hence it suffices to show that 


0 = lim IIQ'Q' T/P'P' -lVr/,zP' 11 = lim \\Q^{Q‘ TfP‘- 1 ‘^Tf, 1 )P^ ||. 
(—>0+ /—>-0+ 

Since and pP are contractions, this follows if 

lim = 0 . 

l-,0+ 


(6.5) 


By the Tietze extension theorem, we may suppose that / is actually defined on M" and has compact 
support there. In particular it will be equicontinuous. Now, to establish (6.5), let g = gi + 92 G 
£^(T/) be arbitrary, where supp gi C T™* and supp 32 C T; \ T™*. By definition, P’'g 2 = 0 so 
QPpfP^g 2 = 0 whereas 

<Zll/||ooArz(T)V2||g2||, 

by the Cauchy-Schwartz inequality. Thus 

\{QpTfP^ - < l‘^\\f\\MPn92\\. 

We now provide estimates for gi. Given x G Si and y G Ti, set 


( 6 . 6 ) 


/(* + y) = ^ 


'\x-x\^<l/2 J\y-y\^<l/2 


f{x + y) dy dx 


and note that 


f{x + y) = j^iQPVfP^dy, 4 


whenever x G and y G . As in the proof of Theorem 6.2 it follows that Q T fP^gi{x) = 
I'^T^ ^gi{x) for x € For such x we thus have 


|(g' TfP^ - fTfMx)\ = \fTj^^^,g,{x)\ < f\\f - /||,.(n,)ll 5 i| 
by Cauchy-Schwartz, and for a; G H \ we get 

liQPPfP^ - l‘^Tf,i)g,{x)\ = ll^^Tf^mix)] < l^^WflUTing.W 


(6.7) 


( 6 . 8 ) 
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due to the definition of QK Combining (6.6)-(6.8) we see that 


< \\{Q^*TfP^ - + |!(Q'*r/P' - < 

- /II0JI51II + Ni{E)^^^l^\\f\Uri\^/^\\g,\\ + |H,|V2z^||/||^iv,(T)i/2||g^||. 


Since S and T are bounded sets, |H/| and |Y/| are bounded by some constant C times and as 
llffill < llsll and II 52 II < II 5 II, it follows that 

||(g'*r/P' - fTf,i)\\ < - /II,2(0,) + Ci/2Ar;(s)i/2/P2||/||^ + ^i/2^d/2||^||^^^(^)i/2^ 

By Proposition |6.3| the last two terms go to 0 as Z goes to 0. The same is true for the first term by 
noting that - /||^ 2 (o,) < ||/ - /||,<x=(o,)Z'^/^| 0 ;|i /2 and 


lim 

i-)- 0 + 


11 / - /|U“(f2,) - 0, 


which is an easy consequence of the equicontinuity of /. 
complete. 


Thereby (6.51 follows and the proof is 


□ 


In particular, we have the following corollary. Note that the domains need not have well-behaved 
boundaries. 

Corollary 6.5. Let T and S be open, hounded and connected domains, and let f be a continuous 
function on cl{Ll). We then have 

Rankrj:= lim RankTyj (6-9) 

(and Rank 0/ = lim;_yo+ Rank @f,i). 

Proof. By Propositions 5.1 and 5.3 in [5, the rank of Ty is independent of T and S. Combining 
this with Theorem |6.2[ it is easy to see that it suffices to verify the corollary for any open connected 
subsets of T and S. We can thus assume that their boundaries are well-behaved. By Theorem |6.4| 
and standard operator theory we have 

Rank Tf < liminf Rank = lim inf Rank Q^Ff /P^ < liminf Rank F, 

1 ^ 0 + 1 ^ 0 + /-)- 0 + 

On the other hand. Theorem |6.2| gives 

lim sup Rank Fy ; < Rank F f. 

i->.o+ 

□ 


7. The multidimensional continuous Caratheodory-Fejer theorem 

In the two final sections we investigate how the PSD-condition affects the theory. This condition 
only makes sense as long as 

S = T, 

which we assume from now on. In this section we show that the natural counterpart of Caratheodory- 
Fejer’s theorem hold for truncated convolution operators 0y on general domains S = T, and in the 
next we consider Fischer’s theorem for truncated correlation operators. 
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Theorem 7.1. Suppose that S = T is open bounded and connected, 12 = S — T, and f G 
Then the operator Qf is PSD and has finite rank K if and only if there exists S and 

Cl,... ,ck > 0 such that 

K 

/ = (7.1) 

k=l 


Proof. Suppose first that Qj is PSD and has finite rank K. By Theorem 4.4 in [3], / is an 
exponential polynomial (i.e. can be written as (6.1 1 ). By uniqueness of analytic continuation, it 
suffices to prove the result for S = T are neighborhoods of some fixed point xq- By a translation, it 
is easy to see that we may assume that xq = 0. Let I assume values 2~fi j G N. For j large enough, 
(beyond J say), the operator has rank K (Corollary 6.5l and Theorem 5.1 applies (upon 


dilation of the grids). We conclude that for j > J the representation (7.1 1 holds (on but the 


^fc’s may depend on j. However, since each grid is a refinement of ^ 2-3 j Proposition |4.1| 

guarantees that this dependence on j may only affect the ordering, not the actual values of the set 
of ^fc’s used in (7.1). We can thus choose the order at each stage so that it does not depend on j. 
Since / is an exponential polynomial, it is continuous, so taking the limit j > 00 easily yields that 
(|!3 holds when a; is a continuous variable as well. 

Conversely, suppose that / is of the form (7.1). Then Qf has rank K by Proposition 4.1 in 
[3] (see also the remarks at the end of Section 2.2). The PSD condition follows by the continuous 
analogue of (5.6). □ 


8. PSD Truncated Correlation Operators 

Theorem 8.1. Suppose that S = T is open bounded and connected, D = S + T, and f G 

The operator P/ is PSD and has finite rank K if and only if there exists G o-nd 

Cl,... ,ck >0 such that 


K 


( 8 . 1 ) 


We remark that the continuous version above differs significantly from the discrete case, even 
in one dimension, since the sequence (A")^^q generate a PSD Hankel matrix for all A G K, whereas 
the base is positive in ( 8 . 1 ). 


Proof. Surprisingly, the proof is rather different than that of Theorem 7.1 First suppose that F/ is 
PSD and has finite rank K. Then / is an exponential polynomial, i.e. has a representation (6.1), by 
Theorem 4.4 in [3]. Suppose that there are non-constant polynomial factors in the representation 
(6.1), say pi. Let N be the maximum degree of all polynomials {pj}j^i. Pick a closed subset S C S 
and r > 0 such that \ S) > 2r. Pick a continuous real valued function g G L^(M'^) with 

support in S that is orthogonal to the monomial exponentials 




(where a G and we use standard multi-index notation), but satisfies ( 5 , = 1 , (that such a 

function exists is standard, see e.g. Proposition 3.1 in 0 ). A short calculation shows that 

- z),gi- -w)} =pi(z-kw)e‘^i'(^+“^ (8.2) 
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whenever |z|,|t(;| < r. Since pi is non-constant, there exists a unit length ly G such that 
q(t) = (ri/t) is a non-constant polynomial in t. Set C = rCi • v- Consider the operator A : 
defined via 

A{(t))= [ (j){t)g{x-rvt). 

^0 

Clearly A*rfA is selfadjoint. It follows by (|8.2[) and Fubini’s theorem that 


{A*rfA{(j)),'tp) = [ f pi{ri't + = 

Jo Jo 

f f q{t + (j){t)dtij}{s)ds. 

Jo Jo 


With h{t) = it follows that the operator F/i : L^([0,1]) —>■ L^([0,1]) is self-adjoint. Either 

by repeating arguments form Section or by standard results from integral operator theory, it is 
easy to see that hit -f s) = h{s + t), i.e. h is real valued. This clearly implies that C G K. Now 
consider the operator B : L^([0,1]) —>■ L^([0,1]) defined by B{g){t) = e~^*g{t). As before we see 
that B*ThB = Vq is PSD. Given 0 < e < 1/2, define C, : L'^{[G,1/‘A) ^ -^-^([0,1]) by C,{g){t) = 
9 (t-f-)-g(t) ^ .^g idgntify functions on [0,1/2] with functions on K that are identically zero 

outside the interval). It is easy to see that 


Cg TgCj — rg-2(q(._|_2j)_2q(.+e)+q(.)), 

which means that also this truncated correlation operator is PSD. Since g is a polynomial, it is easy 
to see that ((7(--|-2e) —2(7(--|-e)-|-(7(-))/e^ converges uniformly on compacts to q”. By simple estimates 
based on the Cauchy-Schwartz inequality (see e.g. Proposition 2.1 in [3]), it then follows that the 
corresponding sequence of operators converges to Tqii (acting on L^([0,1 /2])), which therefore is 
PSD. Continuing in this way, we see that we can assume that q is of degree 1 or 2, where F^ acts on 
an interval [0, 31] where 31 is a power of 1/2. We first assume that the degree is 2, and parameterize 
qit) = a + b{t/l) + c{t/l)‘^. Performing the differentiation trick once more, we see that Fc is PSD on 
some smaller interval, which clearly means that c > 0. Now pick g G L^([0,^]) such that {g, 1) = 1, 
{g, t) = 0, {g, t^) = 0, and consider D : —>■ T^([0, 31]) defined by 


D((co,ci,C 2 )) = cog(-) +cig{- - 1) + C 2 g{- - 21). 


By (8.21, the matrix representation of D*VqD is 


/ 9(0) 

9(0 

9(2/) \ , 

f a 

a + 6 + c 

a + 2b + 4c \ 

M = q(l) 

9(20 

9(3/) = 

a + b + c 

d “h 26 4c 

a -1- 3& -1- 9c j 

V 9(20 

9(30 

9(4/) J ' 

^ Cl H- 26 “h 4c d “h 36 -t- 9c 

a + Ab + 16c J 


which then is PSD. However, a (not so) short calculation shows that the determinant of M equals 
—8c^ which is a contradiction, since it is less than 0 (recall that c > 0). We now consider the case 
of degree 1, i.e. c = 0 and 6^0. As above we deduce that the matrix 


9 (0) 9(0 \^( a a + b + c \ 

q{l) q{2l) J \^a-|-6-|-c a-|-26-|-4c J ’ 


has to be PSD, which contradicts the fact that its determinant is —b^. 
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By this we finally conclude that there can be no polynomial factors in the representation (6.11. 


the form (6.21, i.e. 


f = Et 


By the continuous version of Proposition |4^ (see Proposition 4.1 in [5]), we conclude that / is of 

oCk-X 


for all V G 


conclude that ‘ ^ ^ 
call them henceforth. With this at hand we obviously have 


From here the proof is easy. Repeating the first steps, we 
, by which we conclude that Cfe are real valued. We therefore 


K 






(8.3) 


for all g G F^(S), whereby we conclude that Ck > 0. 

For the converse part of the statement, let / be of the form (8.1). That Tf has rank K has 
already been argued (Proposition 4.1 in |3] ) and that Fy is PSD follows by ( |8.3| ) . The proof is 
complete. □ 


9. Unbounded domains 


For completeness, we formulate the results form the previous two sections for unbounded do¬ 


mains. Fj is defined precisely as before, i.e. via the formula (2.14), except that we now have to 
assume that f{x + -) is in L^(T) for every x G 'E. and vice versa, (see definition 1.1 in [3]). Obvi¬ 
ously, analogous definitions/restrictions apply to 0/ as well. The main difficulty with unbounded 
domains is that exponential polynomials then can give rise to unbounded operators. Following [3], 
we address this by assuming that ft is convex and we let Aid denote the set of directions i) G 
such that the orthogonal projection of O on the half line [0, oo) • i? is a bounded set, and we let 
mt(An) denote its interior. 

Theorem 9.1. Let S = T C be eonvex domains and set = S -|- T. Then Tf is bounded, PSD 
and has finite rank if and only if f is of the form (8.1) and G int{AQ) for all k. 


Proof. This follows by straightforward modifications of the proofs in Section 9 of [3], so we satisfy 
with outlining the details. The “if’ direction is easy so we focus on the “only if”. We restrict 
the operator Ty to functions living on a subset (see Theorem 9.1 [2) to obtain a new operator to 
which Theorem 8.1 above applies. From this we deduce that / locally has the form (8.11. That this 


formula then holds globally is an immediate consequence of uniqueness of real analytic continuation, 
combined with the observation that LI is connected. Finally, the restriction on the ^fe’s is immediate 
by Theorem 9.3 in [3]. □ 


The corresponding situation for truncated convolution operators is quite different. We first note 
that : L^(T) —)■ T^(S) is bounded if and only if Fj : L^(—T) —)■ T^(S) is bounded, as mentioned 
in Section |2.2| and further elaborated on around formula (1.2) in |3]. With this, we immediately 
obtain the following theorem, which was left implicit in [5]. 

Theorem 9.2. Let 2,T C be convex domains and set LI = E — T. Then Qf is bounded and 
has finite rank if and only if f is an exponential polynomial (i.e. fix) = 

Re (j G intiAii) for all j. 


^It was called 0n in [3], see Section 9. 
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However, if now again we let H = T and we additionally impose PSD, the proof of Theorem 


forces 0 = Re Q G mt(An), which can only happen if Aq = K“, since it is a cone. This in turn is 
equivalent to being bounded, so we conclude that 

Theorem 9.3. Let H = T C be convex unbounded domains and set H = S — T. Then Qf is 
bounded and PSD if and only if f = 0. 

10. Conclusions 

Multidimensional versions of the Kronecker, Caratheodory-Fejer and Fischer theorems are dis¬ 
cussed and proven in discrete and continuous settings. The former relates the rank of general 
domain Hankel and Toeplitz type matrices and operators to the number of exponential polynomials 
needed for the corresponding generating functions/symbols. The latter two include the condition 
that the operators be positive semi-definite. The multi-dimensional versions of the Caratheodory- 
Fejer theorem behave as expected, while the multi-dimensional versions of the Kronecker theorem 
generically yield more complicated representations, which are clearer in the continuous setting. Fis¬ 
cher’s theorem also exhibits a simpler structure in the continuous case than in the discrete. We 
also show that the discrete case approximates the continuous, given sufficient sampling. 
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